Li 2 Rh03: A spin-glassy relativistic Mott insulator 
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Motivated by the rich interplays among electronic correlation, spin-orbit coupling (SOC), crystal- 
field splitting, and geometric frustrations in the honeycomb-like lattice, we systematically investi- 
gated the electronic and magnetic properties of Li2Rh03. The material is semiconducting with a 
narrow band gap of A ~78 meV, and its temperature dependence of resistivity conforms to 3D 
variable range hopping mechanism. No long-range magnetic ordering was found down to 0.5 K, 
due to the geometric frustrations. Instead, single atomic spin-glass behavior below the spin-freezing 
temperature (~6 K) was observed and its spin dynamics obeys the universal critical slowing down 
scaling law. First-principles calculation suggested it to be a relativistic Mott insulator mediated by 
both electronic correlation and SOC. With moderate strength of electronic correlation and SOC, 
our results shed new light to the research of Heisenberg-Kitaev model in realistic materials. 

PACS numbers: 71.20.Be, 75.40.Cx, 75.10.Jm, 75.40.Gb, 71.20.-b 
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Ternary transition metal oxides set up a fascinating 
platform for investigating correlated electronic systems. 
Depending on the particular transition metal element 
and crystalline structure, features like high-temperature 
superconductivity (SC) in doped spin-1/2 antiferromag- 
netic (AFM) Mott insulator La2Cu04j 1 - giant negative 
magnetoresistance in La 2 /3Bai/ 3 Mn03 films j 2 - odd par- 
ity SC in Sr2Ru04,— a SC in water intercalated Na x Co02- 
2/H2 0)£ and field induced metamagnetic transition and 
quantum criticality in Sr3Ru 2 0?!i have been observed. 
Electronic correlation is expected to be strongest in 3c? 
transition metals, represented by a small d orbital radius 
and a large Coulomb repulsion U. It weakens as one 
goes from 3d to Ad and 5d transition metals due to the 
spatial extension of d orbits. However, the relativistic 
spin-orbit coupling (SOC) which increases with atomic 
number follows the opposite trend. The recently dis- 
covered exotic non-metal behaviors in those heavy tran- 
sition metal oxides^— remind us of the importance of 
SOC in these materials. One representative example is 
S^IrOzir^ which was confirmed to be a novel Mott insu- 
lator mediated by strong SOC even though the electronic 
correlation is relatively weak, while its structural analog 
Sr 2 Rh04^i shows normal Fermi-liquid metallicity. 

The general formula Li 2 M03 (M = transition metal) 
actually describes two types of crystalline structures: 
Li 2 Mn0 3 -type (C2/m, No. 12)12 and Li 2 Sn0 3 -type 
(C2/c, No. 15) In both crystalline structures, the lay- 
ers of MOq octahedral interstices are alternately filled 
either with Li + only, or with 1/3 Li + and 2/3 M 4+ , 
as depicted in Fig. Q^a), whereas the MOq octahedra 
with M in the center form edge-sharing honeycomb-like 
networkai^[Fig. H{b)]. Tiny difference between these two 
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FIG. 1: (Color online) (a), The crystalline structure of 
Li2Rh03, which is stacked by alternating Li and LiRl^Oe 
layers, (b), Within the LiRl^Oe layer, the RhOe octahedra 
form the honeycomb-like lattice. 



crystalline structures resides in the stacking of Li-M lay- 
ers along c-axis: in the case of Li2Sn03-type, the Sn 4+ 
hexagonal networks in adjacent layers are displaced by 
(0, ± 1/6, 1/2) in lattice coordinates, while in Li2Mn0 3 - 
type, they are displaced by (0, 1/2, 1/2)^ The forma- 
tion of honeycomb-like MOe octahedral network makes 
Li2A/03 a suitable candidate for investigating the in- 
terplays among electronic correlation, SOC, crystal-field 
splitting and geometric frustrations. 

Li2Rh03 and Li2lr03 crystalize in Li2Mn03- and 
Li 2 Sn03-type structures, respectively. Previously, they 
were studied for potential applications as Li-ion battery 
cathode materialsi 14 ' 16 ' 17 More underlying physical prop- 
erties still need to be explored. Herein, we systemati- 
cally studied the electronic and magnetic properties of 
Li 2 Rh03. Our results point out that Li2Rh03 is likely 
to be a spin-glassy Mott insulator with a narrow thermal 
activating gap A ~78 meV, while the spin-freezing tem- 
perature is sample dependent ranging from 5 to 7 K. No 
long range magnetic ordering can be captured in this frus- 
trated system down to 0.5 K. Our experimental results 
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TABLE I: Rietveld refinement of Li2RhC>3. Calculation based 
on space group C2/m (No. 12). The derived a =5.1212(3)1, 
b =8.8469(4)1, c =5.1015(3)1, a= 7 =90°, f3 =109.641(3)°. 
The quality factors of this refinement are R wp =15.40%, 



i? p =11.89%, S =1.79 for SI, and R wp 
S =1.71 for S2. 
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FIG. 2: (Color online) Electrical resistivity of I^RhOs. The 
inset shows resistivity in the Arrhenius plot and 3D-VRH plot. 
The dashed lines are guides to eyes. Thermal activating band 
A ~ 78 meV is estimated in the Arrhenius plot. 



were understood by first principle calculations which con- 
firm the important roles played by both electronic corre- 
lation and SOC. The calculation also pointed out that 
Li2Rh03 is on the boundary of antifcrromagnctically- 
ferromagnetically correlated ground state, which might 
interpret the spin-glassy behavior observed experimen- 
tally. 

Poly-crystalline sample of Li 2 Rh03 was grown by solid 
state reaction method as mentioned elsewhere^ The 
sample quality of Li2Rh03 was checked by X-ray diffrac- 
tion (XRD), performed on a PANalytical X-ray diffrac- 
tometer (Empyrean Series 2) with Cu-K Q i radiation at 
room temperature. Lattice parameters were derived by 
Rietveld refinement on the RIETAN-RF programme.— 
Electrical resistivity and specific heat were measured on a 
Quantum Design physical property measurement system 
(PPMS-9). The dc magnetization measurement was car- 
ried out on a Quantum Design magnetic property mea- 
surement system (MPMS-5) employing both zero-field- 
cooling (ZFC) and field-cooling (FC) protocols. The ac 
magnetic susceptibility was measured on PPMS-9 with 
various frequencies ranging from 10 Hz to 10 kHz. 

The XRD pattern (not shown) guarantees high pu- 
rity of the samples, and all the peaks can be well 
indexed based on the C2/m (No. 12) space group 
iso-structural to Li2MnC>3. The Rietveld refinement 
yields a =5.1212(3)1, b =8.8469(4)1, c =5.1015(3)1, 
a=j =90°, and f3 =109.641(3)°, which are comparable 
with those in the previous literature.— Detailed struc- 
tural parameters are summarized in Table |U The labels 
SI and S2 represent the two Li2Rh03 samples annealed 
at 950 °C and 900 °C, respectively. It should be pointed 
out that there is an antisite disorder between Li + ions 
and Rh 4+ ions. Such anti-site disorder is a common fea- 
ture in Li2-M03 materialS ) 15 i 16 and may has a double-side 



influence on the magnetism: on the one hand, the par- 
tial substitution of M 4+ by non-magnetic Li + ions breaks 
the long range magnetic coupling among M 4+ moments, 
while on the other hand, it also reduces the geometric 
frustrations on the M 4+ honeycomb lattice and thus sta- 
bilizes the magnetic structure. 

The electrical resistivity of SI is shown in Fig. [21 
from which an insulating p vs. T dependence is clearly 
seen. It 300 K, the magnitude of resistivity is 0.48 f2 
cm, about two orders smaller than that of Na2lr03.— 
In the inset of Fig. [3J we show the p vs. 1/T in the 
semi-logarithm plot, and we found that the p{T) curve 
does not well follow the Irrhenius law for thermally- 
activated hopping, viz. p(T) oc cxp(A/T), but is better 
fit to p(T) oc exp[(D/T) 1 / 4 ] which is known as the three 
dimensional variable range hopping (30-VRH) 1 ^ stem- 
ming from the random potential scattering contributed 
by large numbers of defects or disorders. Similar phe- 
nomenon was also observed in Na 2 Ir03^ and we at- 
tributed this to the antisite defect or disorder of Li and 
Rh, as is mentioned in Table HI 1 rough estimate of acti- 
vating energy A ~78 meV will be derived from the data 
200-300 K. The measurement on S2 leads to a similar re- 
sult, lecording to a previous work performed by Todor- 
ova et alyii the p(T) curve obeys the Irrhenius law for 
the temperature region 300-500 K, and the derived en- 
ergy gap is A = 80 meV. This magnitude of A is close 
to our result, and is much smaller than that of Na2lrC>3 
(340 meV)i2£ 111 these confirm that Li2RhC>3 is a narrow 
gap insulator. 

The main frame of Fig. [3]Ja) displays temperature de- 
pendence of dc magnetic susceptibility of SI. For the 
temperatures above 50 K, x{T) obeys Curie- Weiss's law 
[see in the inset of Fig. 03a)], and can be well fitted to 
x(T) = C/(T — 6 W ), where 9w is the Weiss temper- 
ature. The fitting leads to 9w =-50 K. The negative 
6w indicates the dominant 1FM coupling between Rh 4+ 
moments. The fitting also derives the effective moment 
Peff =2.03 \ib (P'B is Bohr's magnon). This magnitude 
of p e ff is close to but relatively larger than the ideal 
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FIG. 3: (Color online) Temperature dependent magnetic susceptibility of Li2Rh03. (a-b), x(T) °f SI and S2 measured under 
various fields, respectively, in both ZFC and FC processes. Inset of (a) displays Curie- Weiss fit of l/x{T) in the high T region, 
while inset of (b) shows hysteretic loop in M(H) below T g . (c), ac susceptibility measurement of SI. Inset of (c): scaling 
plot of log(/) vs. log[(T)/T 9 ) - 1] for SI (solid) and S2 (open), with the best fitted parameters r ~ 10 _u,9(2) s, zv =8.2 (3), 
T g i =5.30 K, and T 9 2 =5.97 K. The solid line is a guide to eyes of this fitting. 
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FIG. 4: (Color online) Main frame: specific heat divided by 
T of Li2Rh03 (red), compared with the lattice contribution 
(black) which was derived by correcting its non-magnetic ref- 
erence Li2Sn03. The solid (open) symbols represent data 
measured under fioH =0 (9 T). Lower inset: C/T of Li2Rh03 
and Li2Sn03 plot in the T 2 scale. Upper inset: magnetic 
contribution to specific heat in Li2Rh03, also shown is the 
magnetic entropy gain S m as a function of T. 

value for low-spin state of Rh 4+ (1.73 /ib) in the case of 
J = S =1/2 and Landc factor g =2, manifesting incom- 
pletely quenched orbital contribution. For the low tem- 
perature region, a sharp peak centered at 5.6 K is clearly 
seen in the curve for /j,qH = 0.01 T. We also observed a 
discrepancy between ZFC and FC modes below this peak 
temperature. With increasing magnetic field, this peak 
looses sharpness and becomes rounded, and meanwhile, 
the discrepancy between ZFC and FC shrinks. This field 
dependent x(^) was confirmed by the isothermal mag- 
netization measurement shown in the inset of Fig. EIb). 



A tiny hysteresis loop with a remanent magnetization 
1.4x10 -3 /zb/R1i and a coercive field 0.6 kOe is evident 
at 2 K. It should be pointed out that these observations 
were reproducible in different batches of samples [e.g., 
the result of S2 is shown in Fig.^b)], although the peak 
position may vary slightly in the range of 5-7 K. 

All these phenomena are hard to be understood by a 
simple AFM or ferromagnetic (FM) transition, but re- 
mind us of the spin-glass transition. We therefore per- 
formed the ac magnetic susceptibility of Li2Rh03, and 
the results are displayed in Fig. [3|c). Indeed, a peak in 
the imaginary part of ac susceptibility x" is evidently 
seen, strongly demonstrating the dissipative process. In 
addition, the real part of ac susceptibility \' shows a peak 
at the freezing temperature (Tf) which shifts towards a 
higher temperature with increasing frequency (/). Such 
frequency dependent dissipative process is a fingerprint 
of a spin-glass transition. The frequency dependence of 
Tf can be described by the conventional " critical slowing 
down" of the spin dynamics^—: 

r(T f ) = T (T f /T g - 1)—, (1) 

where r = 1//, T g is the characteristic temperature of 
spin-glass transition for / — >0, zv is a dynamical expo- 
nent, while To characterizes the intrinsic relaxation time 
of spin dynamics. We show this agreement by plotting 
log(/) vs. log[(T//T g ) - 1] in the inset of Fig. E^c) . The 
results of SI and S2 can be well scaled into a same 
curve in this plot, with the parameters of this scaling 
law zv =8.2(3), tq ~10 _119 ( 2 ) s, while the critical tem- 
peratures for the two samples being T g i =5.30 K and 
T g 2 =5.97 K, respectively. The value of zv is in good 
agreement with the theoretical prediction (7.0-8.0) for 
an Ising spin-glass system. 24 ' 25 The derived tq is in close 
approximation to that of a single atomic spin-glass sys- 
tem which usually possesses a tq in the order of 10" 13 s, 
implying that the observed spin-glass behavior is likely 
to arise from the frustrated single Rh 4+ ions in the hon- 
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FIG. 5: (Color online) The calculated density of state (DOS) and band structure of Li2Rh03, based on (a), LDA; (b), LDA+C/ 
(U =3 eV); (c), LDA+C/+SOC (U =3 eV). The calculations were performed with stripy-AFM structure. 



eycomb lattice, rather than from magnetic domains or 
clusters (for which tq can be as large as 10~ 4 s)i2& 

We now turn to the specific heat of Li2Rh03, as is 
shown in Fig. |4l The measurement was carried out on 
the sample SI. An anomaly is clearly seen at around 7 K. 
Such anomaly differs from a A-shaped specific heat jump 
usually seen in a second order phase transition, indicat- 
ing the absence of long-range magnetic ordering and is 
consistent with the spin-glassy featured Under magnetic 
field, this anomaly is suppressed and a Brillouin-likc po- 
larization trend is observable. Such evolution of specific 
heat under field provides further evidence to the com- 
petition between Zeeman energy and spin-glassy order- 
ing, which is also depicted by the broadened spin-freezing 
peaks in x{T) [Fig. EJa-b)]. We should emphasize that 
no long range magnetic ordering can be captured down 
to 0.5 K in Li2Rh03 by specific heat measurement. For 
comparison, the specific heat of its non-magnetic refer- 
ence Li2SnC>3 was also measured. We fit the specific 
heat of Li 2 Sn0 3 to the formula C Sn /T = y$ n + f3 Sn T 2 , 
and the derived Sommerfeld coefficient is 7^™ =0.18 
mJ/(mol-K 2 ). Such low 7^™ of Li2Sn03 signifies the 
highly insulating electronic property. The slope of this fit 
results in the Debye temperature <dfj l = 418 K. We calcu- 
lated the lattice contribution to specific heat in L12RI1O3 
by correcting Cs n to the molar mass^ and the estimated 
Debye temperature of Li 2 Rh0 3 is eg' 1 = 444 K. The 
magnetic specific heat in Li2Rh03 is obtained by sub- 
tracting the lattice contribution from the total specific 
heat, and the result C m /T is displayed in the upper inset 
of Fig. SI The short range magnetic ordering above T g is 
further represented by the noticeable broad tail in C m /T . 
We then calculated the magnetic entropy gain S m (T) by 
integrating C m /T over T. We found that S m reaches only 
17% of i?ln2 at 20 K, and keeps increasing even for T 
up to 50 K while no evident plateau can be seen (see the 
upper inset of Fig. 2]). There are two sources of entropy 
gain loss. Besides the short range magnetic ordering men- 
tioned above, more magnetic entropy [^70% of i?ln2, 
judging from S m (50K)] should be compensated by the 
residual magnetic entropy stemming from the quantum 



magnetic randomness that persists even at zero temper- 
ature, consistent with the spin-glass scenario. 

To well understand these experimental results, we per- 
formed first-principles calculation^ For the id transi- 
tion metal element, both electronic correlation and SOC 
should be taken into account. The calculated den- 
sity of state (DOS) and band structure are shown in 
Fig. 03 We started with the local density approximation 
(LDA), from which we derived a metallic electronic state 
[Fig. 03a)] as Rh 4+ has the half filled ionic configuration 
of 4d 5 . The application of Coulomb repulsion U =3 eV— 
greatly reduces the DOS at Fermi level, however, there 
are still two bands crossing the Fermi level [Fig. [5fb)] 
and forming a semi-metal-like band structure. We should 
point out that such semi-mctal-like band structure is ro- 
bust to Coulomb repulsion and will persist even under 
U >4 eV (data not shown). SOC was then employed, 
and the combination of U and SOC successfully elimi- 
nates the band crossing at Fermi level and thus opens 
a gap A ~65 meV in the DOS spectrum, as shown in 
Fig. 03c). This magnitude of the energy gap is close to 
the thermal activating gap (78 meV) derived experimen- 
tally. The calculated SOC splitting is ~10 meV, much 
smaller than that of Na2lr0 3 ^ and the orbital- and spin- 
moment contributions are respectively (L) =0.21 /^/Rh 
and (5)=0.19 ^s/Rhi^L We should also point out that 
merely SOC can not open a gap at Fermi level, either 
(data not shown). Therefore, Li2Rh0 3 is suggestively to 
be a relativistic Mott insulator driven by both electronic 
correlation and SOC. 

The calculation also helps us understand the spin- 
glassy feature of Li2Rh03. The study of such SOC me- 
diated honeycomb lattice turns to the Hcisenbcrg-Kitaev 
modclj22r— in which FM ordering, AFM ordering with 
Neel-/ stripy-/ zigzag-type, or spin-liquid state emerges 
depending on the particular anisotropic magnetic cou- 
plings. In the case of Li2Rh03, the calculation based on 
LDA+U prefers a FM ground state, but is challenged 
by several other magnetic configurations. When SOC is 
turned on, the ground state switches to stripy- or zigzag- 
type AFM configuration, and still many other magnetic 
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configurations are comparable in energy^ In this situa- 
tion, perturbations such as disorders or defects (see Table 
|l| are likely to change the magnetic ground state. It is the 
fact that Li2Rh03 embeds in a regime close to the multi- 
phase boundary that results in the spin-glass nature. In 
addition, according to Choi et al's result of inelastic neu- 
tron scattering experiment on single crystalline Na2lr03, 
there is some proportion of stacking fault of well-ordered 
honeycomb layers along c axis ; 37 ' 38 which is hardly re- 
solvable by powder XRD pattern^ Such stacking fault 
might also appear in Li 2 Rh03 and account for the spin- 
glass ordering. To clarify the magnetism of Li2Rh03, sin- 
gle crystals are highly needed. With moderate strength 
of electronic correlation and SOC, our result sheds new 
light to the research of Heisenberg-Kitaev model in re- 
alistic materials and calls for more investigations in the 
future. 

To summarize, we systematically studied the electronic 



and magnetic properties of Li2RhC>3 on poly-crystalline 
samples. Our experiment confirms that Li 2 Rh03 is a 
spin-glassy insulator with a narrow gap A ~78 meV. 
This picture is supported by first-principles calculation 
which verifies the combination of electronic correlation 
and SOC. The calculation also points to many nearly 
degenerated magnetic configurations, which possibly il- 
lustrates the spin-glass behavior. Our result provides a 
unique case for the studies of Heisenberg-Kitaev model 
in realistic materials with moderate electronic correlation 
and SOC. 

Y. Luo thanks Hui Xing and Ying Liu for helpful dis- 
cussions. This work was supported by the National Basic 
Research Program of China (Grant Nos. 2011CBA00103, 
2012CB927404 and 2010CB923003) , the National Science 
Foundation of China, and the Fundamental Research 
Funds for the Central Universities of China. 



1 J. G. Bednorz, and K. A. Miiller, Z. Physik B64, 189 
(1986). 

2 R. von Helmolt, J. Wecker, B. Holzapfel, L. Schultz, and 
K. Samwer, Phys. Rev. Lett. 71, 2331 (1993). 

3 Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, T. Fu- 
jita, J. G. Bednorz, and F. Lichtenberg, Nature 372, 532 

(1994) . 

4 K. D. Nelson, Z. Q. Mao, Y. Maeno, and Y. Liu, Science 
306, 1151 (2004). 

5 K. Takada, H. Sakurai, E. T. Muromachi, F. Izumi, R. A. 
Dilanian, and T. Sasaki, Nature 422, 53 (2003). 

6 R. S. Perry, L. M. Galvin, S. A. Grigera, L. Capogna, A. 
J. Schofield, A. P. Mackenzie, M. Chiao, S. R. Julian, S. 
I. Ikeda, S. Nakatsuji, Y. Maeno, and C. Pfieiderer, Phys. 
Rev. Lett. 86, 2661 (2001). 

7 D. Yanagishima, and Y. Maeno, J. Phys. Soc. Jpn. 70, 
2880 (2001). 

8 A. S. Erickson, S. Misra, G. J. Miller, R. R. Gupta, Z. 
Schlesinger, W. A. Harrison, J. M. Kim, and I. R. Fisher, 
Phys. Rev. Lett. 99, 016404 (2007). 

9 Y. Singh, and P. Gegenwart, Phys. Rev. B 82, 064412 
(2010). 

10 B. J. Kim, Hosub Jin, S. J. Moon, J. Y. Kim, B. G. Park, 
C. S. Leem, J. Yu, T. W. Noh, C. Kim, S. J. Oh, J. H. 
Park, V. Durairaj, G. Cao, and E. Rotenberg, Phys. Rev. 
Lett. 101, 076402 (2008). 

11 R. S. Perry, F. Baumberger, L. Balicas, N. Kikugawa, N. 
J. C. Ingle, A. Rost, J. F. Mercure, Y. Maeno, Z. X. Shen, 
and A. P. Mackenzie, New J. Phys. 8, 175 (2006). 

12 M. Jansen, R. Hoppe, Z. Anorg. Allg. Chem. 397, 279 
(1973). 

13 G. Lang, Z. Anorg. Allg. Chem. 348, 246 (1966). 

14 V. Todorova, and M. Jansen, Z. Anorg. Allg. Chem. 637, 
37 (2011). 

15 H. Kobayashi, R. Kanno , Y. Kawamoto, M. Tabuchi, 
O. Nakamura, and M. Takano, Solid State Ionics 82, 25 

(1995) . 

16 H. Kobayashi, M. Tabuchi, M. Shikano, H. Kageyamaa, 
and R. Kannob, J. Mater. Chem. 13, 957(2003). 

17 M. J. O'Malley, H Verweij, and P. M. Woodward, J. Solid 



State Chem. 181, 1803 (2008). 

18 F. Izumi, and K. Momma, Solid State Phenom. 130, 15 
(2007). 

19 P. V. E. McLintock, D. J. Meredith, and J. K. Wigmore, 
Matter at low temperatures, (Blcakie, Glasgow, 1984), pp. 
82. 

20 R. Comin, G. Levy, B. Ludbrook, Z. H. Zhu, C. N. Veen- 
stra, J. A. Rosen, Yogesh Singh, P. Gegenwart, D. Strieker, 
J. N. Hancock, D. van der Marel, I. S. Elfimov, and A. 
Damascelli, Phys. Rev. Lett. 109, 266406 (2012). 

21 K. Binder, and A. P. Young, Phys. Rev. B 29, 2864 (1984). 

22 J. A. Mydosh, Spin Glasses: An Experiment Introduction 
(Taylor & Francis, London, 1993). 

23 K. Gunnarsson, P. Svedlindh, P. Nordblad, L. Lundgren, 
H. Aruga, and A. Ito, Phys.Rev. Lett. 61, 754 (1988). 

24 S. Kirkpatrick, and D. Sherrington, Phys. Rev. B 17, 4384 
(1978). 

25 I. A. Campbell, Phys. Rev. B 33, 3587 (1986). 

26 A. K. Bera, and S. M. Yusuf, Phys. Rev. B 86, 024408 
(2012). 

27 K. H. Fischer, and J. A. Hertz, Spin Glasses (Cambridge 
University Press, Cambridge, 1991) p. 301. 

28 M. Bouvier, P. Lethuillier, and D. Schmitt, Phys. Rev. B 
43, 13137 (1991). 

29 The calculations were performed with plane-wave basis 
projected-augmented wave code VASP, with the PBE fla- 
vor of GGA chosen to be the exchange-correlation func- 
tional. Total energy cutoff was set to be 540 eV and 8x5x8 
T-centered K-mesh was used. For more details of the cal- 
culation, please refer to C. Cao, Y. Luo, Z. Xu, and J. Dai, 
arXiv: 1303.4675 (2013). 

30 S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. 
Humphreys and A. P. Sutton, Phys. Rev. B 57, 1505 
(1998). 

31 X. Wan, A. M. Turner, A. Vishwanath, and S. Y. Savrasov, 
Phys. Rev. B, 83, 205101 (2011). 

32 A. Kitaev, Ann. Phys. (NY.) 321, 2 (2006). 

33 J. Chaloupka, G. Jackeli, and G. Khaliullin, Phys. Rev. 
Lett. 105, 027204 (2010). 

34 H. C. Jiang, Z. C. Gu, X. L. Qi, and S. Trebst, Phys. Rev. 



B 83, 245104 (2011). 

X. Liu, T. Berlijn, W. G. Yin, W. Ku, A. Tsvelik, Y. J. 
Kim, H. Gretarsson, Y. Singh, P. Gegenwart, and J. P. 
Hill, Phys. Rev. B 83, 220403(R) (2011). 
Y. Singh, S. Manni, J. Reuther, T. Berlijn, R. Thomale, 
W. Ku, S. Trebst, and P. Gegenwart, Phys. Rev. Lett. 108, 
127203 (2012). 

S. K. Choi, R. Coldea, A. N. Kolmogorov, T. Lancaster, 
I. I. Mazin, S. J. Blundell, P. G. Radaelli, Y. Singh, P. 
Gegenwart, K. R. Choi, S. W. Cheong, P. J. Baker, C. 
Stock, and J. Taylor, Phys. Rev. Lett. 108, 127204 (2012). 
F. Ye, S. Chi, H. Cao, B. C. Chakoumakos, J. A. 
Fernandez-Baca, R, Custelcean, T. F. Qi, O. B. Korneta, 
and G. Cao, Phys. Rev. B 85, 180403(R) (2012). 



